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Abstract 

We show that faint red stars do not contribute significantly to the mass budget 
of the Galaxy or to microlensing statistics. Our results are obtained by analyzing 
two long exposures of a high-latitude field taken with the Wide Field Camera on 
the newly repaired Hubble Space Telescope. Stars are easily distinguished from 
galaxies essentially to the limiting magnitudes of the images. We find five stars 
with 2.0 < V - I < 3.0 and / < 25.3 and no stars with V - I > 3.0. Therefore, 
main-sequence stars with Mj > 10 that are above the hydrogen-burning limit 
in the dark halo or the spheroid contribute < 6% of the unseen matter. Faint 
red disk stars, M-dwarfs, contribute at most 15% to the mass of the disk. We 
parameterize the faint end of the cumulative distribution of stars, $, as a function 
of luminosity Ly, / din Ly oc Ly'^ . For spheroid stars, 7 < 0.32 over the range 
6 < My < 17, with 98% confidence. The disk luminosity function falls, 7 < 0, for 
15 ^ My ^ 19. Faint red stars in the disk or thick disk, and stars with My < 16 
in the spheroid contribute r < 10 ~® to the optical depth to microlensing toward 
the Large Magellanic Cloud. 
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1. Introduction 

Historically, the direct detection of faint stars has been limited by the diffi- 
culty in distinguishing with ground-based observations faint images of stars from 
galaxies, which are much more numerous at faint magnitudes (V ^ 20). Here we 
present the results of a search for faint red Galactic stars, as originally suggested by 
Bahcall, Guhathakurta, & Schneider (1990), in deep Wide Field Camera (WFC2) 
images taken with the newly repaired Hubble Space Telescope (HST). The images, 
which greatly improve previous limits on the number of faint red stars, cover 4.4 
square arcmin and total 2.00 hr in V' (F606W) and 2.83 hr in /' (F814W). The 
field contains 5 stars with 2.0 <V — I < 3.0 and no stars with V — I > 3.0, down 
to the the magnitude limit / < 25.3. The field analyzed in this Letter = 241, 
h = —51), has the longest exposures among seven HST WFC2 fields that we are 
currently investigating. In a future paper, we will report our techniques for ana- 
lyzing the other fields and the resulting measurements of the luminosity functions 
(LFs) of the disk, thick disk, and spheroid. However, we have checked one of the 
other fields = 148, b = —57), which has a magnitude limit / < 24.7, in order 
to make certain that the primary field is not anomalously devoid of stars. This 
second field contains five stars with 2.0 < V — I < 3.0 and none with V — I > 3.0. 
It is therefore consistent with the primary field. 

The resolution of these images is set by the pixel scale, 0."l. One can distin- 
guish between stellar and extended objects down to the magnitude limit, several 
magnitudes deeper than all previous star-count studies. (See Bahcall 1986 for a 
review and e.g., Majewsky 1992; Richer & Fahlman 1992; Richstone et al. 1992, 
Gould, Bahcall, & Maoz 1993 for several deeper and more recent studies using a 
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variety of techniques.) For faint red stars, we establish improved constraints on 
contributions to the mass of the Galaxy, the slopes of LFs, and the optical depths 
to microlensing. The limits on faint red stars obtained from the HST WFC2 star 
counts are stronger than either the dynamical or microlensing limits for all four 
Galactic components considered in this paper. We cannot put significant limits on 
brown dwarfs because they are too faint. 

For purposes of this Letter, we define the faint endpoint of the main-sequence 
as the reddest known nearby star in V — I, i.e., Mj = 14, V — I = 4.7 (see 
Monet et al. 1992). (Fainter stellar objects are known, but these are also bluer.) 
We refer somewhat loosely to this color-dependent endpoint as the "hydrogen- 
burning limit". The actual magnitude of the faintest hydrogen-burning stars is 
uncertain, but is probably close to the value adopted here. For the disk, we also 
give results in terms of a hydrogen-burning limit of Mj = 15 (A. Burrows 1994, 
private communication). 

2. Observations 

We use 6 V and 7 /' WFC2 CCD iJ^T observations of a field at / = 241.1237, 
b = —51.2286 from the Guaranteed Time Observer's parallel observing program. 
These exposures ranged from 500 s to 2100 s, for a total of 7200 s in V' and 10,200 
s in /'. The individual images are aligned to better than 0.1 pixel. We combined 
the images by calculating their weighted and scaled means, but for each pixel we 
excluded in sequence the highest a outliers until each included pixel lay within 3 a 
of the mean of the others. This removed the cosmic-ray events efficiently. We then 
removed the dead pixels and hot pixels (as determined by median stacking images 
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of 10 other fields) by a quadratic fit to the eight nearest pixels. This procedure 
produced combined images with background noise in electrons per second of 3.3 X 
10~^ e pixel"-^ in /' and 4.7x10"^ e s~-^ pixel"-^ in V' . We constructed a catalog 
of stars in two stages. First, we established a list of stellar candidates on each of 
the resulting (800 X 800) WFC2 images in each of the two bands by searching (by 
eye) for objects that are not obviously extended and that satisfy one or more of the 
following criteria: I) V - I > 2.0 and /' flux > 0.25 e s'^ 2) /' flux > 0.39 e s'^ 
or 3) V' flux > 0.54 e . Second, we subjected each candidate to a star/galaxy 
discrimination test. To construct this test, we examined an image of a fleld at low 
galactic latitude (6 = 7°) with > 20 bright stars on each chip. We compared the 
radial proflles of all the stars from each chip, which showed that the point spread 
function is nearly constant over each chip. We used this superposition of radial 
proflles as a test template for deciding if objects are stars. All of the candidates 
fall clearly into 'stellar' and 'non-stellar' categories. Finally, we repeated the entire 
procedure replacing the 'by-eye' search for candidates with an automated search. 
The resulting catalog of stars was identical to the 'by-eye' search. 

Any stellar object detected in l' had a measurable flux in V' and vice- versa. 
The main results of this Letter concern red stars; hence, the important flux limit 
is in /'. In terms of Johnson/Cousins /, the limiting magnitude is / < 25.3 for 
2 < V — / < 5. Although it is not important for this Letter, the limiting V 
magnitude is V ^ 25.6 + 0.3(V — /) for < V — / < 5. The statistical errors in the 
photometric measurements and the (usually-larger) flat-fleld uncertainties (^ 0.1 
mag) are too small to affect our conclusions. The list of stellar candidates may 
include faint compact galaxies as well as stars. However, we are confldent that we 
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have not excluded any stars. 

We transformed the V' , I' system to standard Johnson/Cousins V and / band- 
passes as follows. We convolved the spectrum of each of 175 Gunn & Stryker 
(1983) standard stars with V and / throughput curves obtained from D. Schneider 
(private communication 1994) and with the V' and /' combined filter and chip 
response curves (Burrows et al. 1993). The V' and /' response curves, which differ 
slightly from chip to chip, were made available to us in electronic form by the 
Space Telescope Science Institute. We established the V and / zero points using 
the preliminary ground-based calibration by J. Holtzman (private communication 
1994). These zero points agree to within ^0.1 mag with the values expected based 
on the HSTWFC2 manual (Burrows et al. 1993). We fit the Gunn & Stryker stars 
to a piece-wise linear relation to obtain, 

Vn =OLV - 2.5 \og{Cv') + ISviy + Sn, 

(2.1) 

4 =ai - 2.51og(C70 + f3iiV + Sn, 
for standard Johnson/Cousins V and / on chips n = 2, 3, 4 as functions of the Cyi 
and Cji, the fluxes in electrons per second. We find ay = 24.90, f3v = 0.284, 
aj = 23.79, /3/ = -0.033 for < V - / < 1.90; ay = 25.20, = 0.126, 
aj = 23.72, /3/ = 0.004 for 1.90 < V - I < 3.90; ay = 24.14, ^y = 0.400, 
aj = 23.17, /3/ = 0.142 for 3.90 < V - I < 5.0; and S2 = -0.03, ^3 = 0.09, 
^4 = -0.06. 

Figure 1 shows the color-magnitude diagram of stellar objects in the field; the 
reddest object has V — I = 2.9. 

The star counts for the disk and thick disk were predicted using the color- 
magnitude relation Mj = 3.45 + 2.25{V - I) for 2.0 < - / < 4.7. For spheroid 
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and halo stars, we use Mj = 7.25 + 1A{V - I) for 2.0 <V-I < 4.5. The disk and 
thick relation is obtained from data for dwarf stars given by Monet et al. (1992). 
Conservatively, the color-magnitude diagram for the spheroid and halo stars was 
estimated from the lower envelope of the subdwarfs. 

3. Comparison With Dynamical Limits On Dark Matter 

If the dark halo of the Galaxy were made of main-sequence stars Mj = 14, 
they would be visible in the iJST field to a distance d = 1.8 kpc. We would expect 
iV = 65 such stars in the field (assuming a local halo density of 9 X lO~^iW^0pc~^ 
Bahcall et al. 1983). No stars are observed with V — I > 3.0 and only 5 stars 
(~ 4 of which are probably in the disk or thick disk based upon their inferred 
distance moduli) are detected with 2 < V — I < 3. If there are halo stars at the 
main-sequence limit, Mj = 14, they constitute < 6% of the dark halo with 98% 
confidence. Making the very conservative assumption that all of the observed stars 
with V > 2 are in the halo, we find that stars in the range 10 < Mj < 13 contribute 
< 5% of the dark halo and stars with 10 < Mj < 14 contribute < 6%. 

The same limits on the fractional contributions to the mass of the Galaxy as 
were given above for the halo also apply to spheroid stars. 

The maximum local column density of mass in the disk that is consistent with 
the observed rotation curve is Sm^x ~ 

135M0PC-2, while the observed density of 
visible matter. Sobs = 48M^0pc~^ (Bahcall 1984). Various studies of the vertical 
motions of tracer stars further constrain the 2 a upper limit on the total mass (disk 
+ halo) within ~ 1 kpc of the plane to TOiW^Q pc~^ to 85iW^0pc~^, depending on 

6 



the technique used (Kuijken & Gilmore 1991; Bahcall, Flynn, & Gould 1992; and 
Flynn & Fuchs 1994). Here, we place independent limits on the contribution to the 
total column density by faint red stars, of mass m and absolute magnitude Mj, in 
a thin disk or a thick disk. 

For an exponential disk of scale height h, the total number of stars counted, 
iV, per unit area, $7, toward Galactic latitude 6, is 



d 

dN S / , 9 / rsin6\ ^h^ csc^ b J d sin b\ , , 



dVl 2hm 





where ((x) = 1 - (1 + 3; + /2) exp(-3;) and d = dex(0.2[/iiin - Mj] + 1) pc. Note 
that C(3) ~ 0.6 and C(^) ~ 1 for x ^ 5 (luminous stars). For an ensemble of stars, 
we may replace the quantity ^/m by its average, ^/m {(/m). 

For V — I > 3.0, we detect iVg^s = stars. We can therefore rule out iV > 4 
at the 98% confidence level and hence place an upper limit of 



S < 8M0pc-2( ( O.2M0'J . (3.2) 



This corresponds to a limit of S < 8 iW^© pc~^ for stars with 10.2 < Mj < 13.5 
where ^ ~ 1. This would correspond to < 15% of the total mass of the disk. If 
all the stars have Mj = 14, the limit actually becomes stronger since the stars 
have lower mass. Taking account of this and the lower we find S < 5 il^© pc~^ 
(for Mj = 14) and S < 8M0PC-2 (for Ml = 15). Previous estimates gave 
S ~ I6M0PC-2 for all thin disk stars with My > 6.8, i.e., - / > 1.1 (Bahcall 
1984). 



For a thick disk of scale height h ~ 1.4 kpc (Gilmore & Reid 1983), we obtain 
98% confidence Hmits of S < 3.3M0pc-2 for stars with - / > 3.0 (M/ > 10.2) 
and brighter than Mj < 14 and S < IOMqpc"^ for Mj < 15. 



4. Constraints On Disk and Spheroid Luminosity Functions 

The faint end of the disk and spheroid luminosity functions (LFs) have in the 
past been measured primarily from nearby stars. Important information about 
the LFs has been obtained from star counts down to V ~ 20, beyond which 
star/galaxy confusion makes further progress difficult with ground-based images. 
However, nearby samples suffer from several limitations. For the spheroid, the 
absolute number of faint nearby stars is small. Moreover, kinematically selected 
spheroid samples may be biased. For the disk, a much larger fraction of nearby 
stars are young than is characteristic of the disk population as a whole. For the 
faintest stars, the volume probed by the HST WFP2 field is equivalent to a Aiv all 
sky sample complete to / = 13, much deeper than any local volume survey. Star 
counts from this field therefore provide the best constraints yet on the shape of the 
faint end of the LFs. 

We establish constraints by comparing the predicted number of stars, iV, in the 
region V — I > 3.0, / < 25.3 with the observed counts, iVg^s = 0. We parameterize 
the faint end by a power-law exponent, 7, with the cumulative number of stars, $, 
growing as a function of V-band luminosity, / din Ly oc Ly'^ . 

For the disk, we find that a conventional Bahcall-Soneira model (Bahcall 1986) 
with a fiat (7 = 0) LF beyond My = 15 (where the LF was previously poorly 
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constrained) predicts 4.9 stars and is thus ruled out with 99% confidence. The 
disk LF therefore falls for 15 > My > 18.7, i.e., to the hydrogen burning limit. 

For the spheroid, the LF has been poorly constrained beyond My = 6. Nor- 
malizing the spheroid LF locally to 1/500 (Bahcall, Schmidt, and Soneira 1983) of 
the disk LF at My = 6, we find that power laws with 7 > 0.32 are ruled out at the 
98% confidence level. If the spheroid is normalized locally at 1/800 of the disk LF 
(as suggested by Gilmore & Reid 1983), then 7 > 0.37 is similarly ruled out. The 
HST ^e\d effectively probes the spheroid LF over a mag range 14 ^ My ^ 17. 

The upper limits obtained here disagree with the average value, 7 = 0.63, found 
by Richer & Fahlman (1992) from faint star counts over the interval 6 < My < 14. 
Richer & Fahlman argue that their result showed that the spheroid LF is consistent 
with that found in some globular clusters. However, using their value of 7 = 0.63 
and a local normalization of 1/500 at My = 6, we calculate we should have detected 
iV ~ 75 stars with V — I > 3.0. in the iJST field, against none observed. We should 
also have detected iV ~ 45 with 2<y-/<3(12< My < 14.5) against iVobs < 5 
observed. The large slope obtained by Richer & Fahlman may result from the fact 
that their normalization of the spheroid at My = 5-6 is 10~® stars mag~'^ pc~^, two 
orders of magnitude below the measured value (see Bahcall, Schmidt, & Soneira 
1983). Also, if the local spheroid denisty is as large as proposed by Caldwell & 
Ostriker (1981), 1.1 x 10~^ Mqpc~^ , then this mass must be primarily in objects 
fainter than My ~ 16. 
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5. Microlensing and Star Counts 

Two groups have recently detected a total of 6 candidate microlensing events 
toward the LMC. (Alcock et al. 1993, MACHO; Alcock et al. 1994; Aubourg et 
al. 1993, EROS). These events were found, as originally suggested by Paczyhski 
(1986), by monitoring more than 10 million LMC stars over about 300 days, in 
order to detect massive compact halo objects (Machos) that may make up the 
dark matter in galactic halos. 

On the basis of the 4 detections by MACHO, we estimate a minimum optical 
depth of Tmiii = 4 X 10~® for events long enough to have been observed a number 
of times (and assuming a / = 100% detection efficiency and that all candidates are 
true microlenses). The true optical depth is probably higher because / < 100%. 
Several workers have suggested that the events are due to a faint, previously unseen, 
component of the disk, thick disk, or spheroid (Gould, Miralda-Escude, & Bahcall 
1994; Gould 1994a; Giudice, MoUerach, & Roulet 1994). Here we show that low- 
mass stars (late M dwarfs or subdwarfs) do not contribute significantly to the 
observed microlensing events. 

Consider a population of stars that are distributed exponentially in height 
above the Galactic plane with scale height h and with total column density S. 
From equation (3.1) and Gould (1994a), we find the number of star counts in an 
$7 = 4.4 square arcmin field required to account for a given the optical depth, r, 
toward the LMC is 

2tvG \smH\ V350pcyV\"^/ 7 V^miny 

where for the sake of definiteness we have written the result in terms of t^[j^ = 
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4 X 10 . Because the LMC and the HST field are both at relatively high latitude, 
we have ignored the radial dependence of the disk density. 

Equation (5.1), combined with the fact that we detect no stars with V — I > 3.0, 
implies that thin-disk red stars Mj < 13.5, bright enough to be seen at a distance 
corresponding to five disk scale heights d = 5hcscb = 2.3 kpc are responsible for 
< 0.08 

'''min) with 98% confidence. If all the stars have Mj — 14, the maximum 
contribution is smaller, < 0.05 Tmin (< 0.08 Tmin at Ml = 15). 

For a thick disk of scale height h ~ 1.4 kpc (as suggested by Gilmore & Reid 
1983), red stars with Mj < 14 contribute < 0.12rinin- The limit rises to < 0.36rinin 
for Ml =15. 

Gould (1994b) has calculated the optical depth to a source lying within a 
spheroidal mass distribution with radial profile as seen by an infinitely distant 
observer. These same formulae apply for an observer within a spheroidal mass 
distribution and looking at an infinitely distant source. Since the LMC is at a 
finite distance, there is a ~ 20% reduction in r. 

The number of stars of mass m and absolute magnitude Mj observed per unit 
area toward Galactic coordinates / and b is 

d{Mi) 

dN = P^ I dr (52) 

dO. m J (r2 + i22 _ 2ri2o cos cos 6)^2 ^ ' 



The expected number of stars in the HST field per unit optical depth toward 
the LMC rises steeply with the index, n. We investigate the n = 3 model in order 
to place limits on all steeper spheroid models. For a poplulation of spheroid stars 
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with V - I = 3.0 {Mi = 11.5) we expect to detect N = 70t/t^[^ stars versus 
iVobs < 1 detected with V -I > 2.6. At - / = 2.5 and - / = 2.0, the number 
of expected stars rises to N = 14:0t / t^[j^ and N = 270t /t^[j^ respectively, while 
the total number of stars detected rises to iVg^s < 5. These calculations imply that 
spheroid stars 10 < Mj < 12 account for < O-OdT^in with 98% confidence. Stars 
with 12 < Mj < 13.5 account for < l-Or^in at the 98% confidence level. 

6. Summary and Discussion 

Faint red stars down to the hydrogen- burning limit do not contribute signifi- 
cantly to the mass budget of the Milky Way. Faint red stars in the dark halo or 
in the spheroid constitute < 6% of the mass of the Galaxy. M-dwarfs in the disk 
constitute less than < 15% of the disk mass. 

Late M dwarfs in the disk or in a thick disk also do not contribute significantly 
to the microlensing events that have been detected toward the LMC. If main se- 
quence stars in the spheroid are responsible for these microlensing events, they 
must be fainter than Mj = 13. 
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FIGURE CAPTIONS 



olor-magnitude diagram of stellar objects detected in a 4.4 square arcmin field 
toward / = 241 b = —51, using the repaired HST WFC2. Standard John- 
son/Cousins V and / magnitudes were obtained from the instrumental mag- 
nitudes using equation (2.1). The mean magnitude /-band limit for the 3 
chips (solid line) was evaluated by subsituting the flux limits (see text) into 
equation (2.1) with Sn = 0. Objects above the flux limit in either V ot I are 
included in the study. 
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